The thermal conductance of a two-channel-Kondo model 
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A theory of thermal transport in a two-channel Kondo system, such as the one formed by a small 
quantum dot coupled to two leads and to a larger dot, is formulated. The interplay of the two 
screening constants allows an exploration of the Fermi liquid and non-Fermi liquid regimes. By 
using analytical, as well as numerical renormalization group methods, we study the temperature de- 
pendence of the thermal conductance and the Lorentz number. We find that in the low-temperature 
limit, the Lorentz number attains its universal value, irrespective of the nature of the ground state. 

PACS numbers: 72.10.Fk, 73.63.Kv 



I. INTRODUCTION 

Quantum dots (QD) have long been recognized by ex- 
perimentalists and theorists alike for their ability to em- 
body several important theoretical paradigms. Some of 
the most important refer to the possibility to replicate 
in a tunable structure the Kondo effect. First associ- 
ated with the anomalous resistivity values in metals, the 
Kondo effect originates in the antiferromagnetic interac- 
tion of conduction electrons with a local magnetic mo- 
ment of spin 1/2. In this respect, the degenerate ground 
state of the electron liquid in the quantum dot functions 
as a magnetic impurity, while the electrons in the leads 
connected to the dot behave as the surrounding normal 
metal^. The multi-channel Kondo problem involves two 
or more electron modes in the screening of the magnetic 
impurity. The QD representation of this situation has 
been realized several years ago and is based on the dif- 
ferent interactions established between an electron in a 
small dot, that functions as a magnetic impurity, with 
the electrons in two leads and with those located in a 
larger doli 2 -. The relative strength of the screening real- 
ized by these two channels determines the characteristic 
behavior of the system. 

If in the case of a single-channel Kondo (1CK) effect, 
at temperatures below the characteristic Kondo tempera- 
ture, Tk, the system behaves like a Fermi liquid (FL) and 
is characterized by analytic functions of temperature, the 
interplay between the two screenings in the case of two 
channel Kondo (2CK) problem, permits the realization of 
a FL, when one screening dominates the other one, or of 
a non- Fermi liquid (NFL), when the screenings are equal. 
Then each reservoir is trying to screen the magnetic im- 
purity, but in the symmetric limit J\ = Ji-, no one suc- 
ceed. This leads to a NFL type of ground state, where, in 
contrast to the single-channlc Kondo problem, the local 
impurity is only partially screened. The impact of the 
Kondo correlations on electronic transport has long been 
an area of active research. Considerably less is known 
about their effect on thermal transport, which only re- 
cently has been discussed in several theoretical^— and 
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FIG. 1: (color online) The universal behavior of the thermal 
conductance of a 2CK model (corresponding to channel 1). 
Three regimes are clearly visible when J\ 7^ J2. For T <C 
T* the system is in the Fermi liquid regime, followed by a 
crossover (T* < T < T K ) to the perturbative (T > T K ) 
regime. The non- Fermi fixed point is realized for J\ = 3i. a 
and 7 are two numerical constants of order 1, and Go is the 
universal conductance Go = 2e 2 /h. 



experimental^ studies that focus on the 1CK model. In 
this paper we extend such considerations to the study of 
the thermal conductivity of a double quantum dot sys- 
tem that can exhibit a two-channel Kondo stated. In 
this picture, the single electron in the small dot interacts 
with the electrons in the leads, with a screening integral 
Ji and with the electrons in the bigger dot with a screen- 
ing integral J2. Besides the Kondo temperature itself, 
Tk, the relationship between J\ and J2, expressed by 
K = {J\ — J2)/ J 2 sets another energy scale of the prob- 
lem through the characteristic temperature T* = K 2 Tk, 
smaller than Tk- For the particular experimental setup 
realized in Ref. [13 the Kondo temperature ranges be- 
tween 30 and 130 mK. 

In an approach that combines analytical arguments 
with Wilson's numerical renormalization group (NRG) 
method, we calculate the thermal conductivity and the 
Lorentz number of this system, and investigate their be- 
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havior as a function of temperature. In this algorithm, 
the thermal conductance k(T) corresponding to channel 
1 is obtained as 
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where 6*(A') is the Hcaviside function and sgn(fc) is the 
signum function; a and 7 are two numerical constants^ 
of order 1. A sketch of this result is presented in Fig. [T] 
We find that when Jj ^ J 2 , at temperature T < T* , the 
ground state of the system is a Fermi liquid (FL) and 
the thermal conductance varies as k(T)/T oc (T/T*) 2 . 
The symmetric channel problem, when J\ = J2, corre- 
sponds to a NFL behavior when the only energy scale in 
the problem is set by T K . Then, k{T)/T oc y/T/T K . 
In the large temperature regime, when T Tk, a 
1/log 2 (T/T K ) behavior is obtained for k(T)/T. 

The nature of the ground state strongly affects the 
Lorentz number, namely, L(T) = k(T)/TG(T), with 
G(T) the charge conductance. In the Fermi liquid regime, 
quasiparticles transport both charge and energy and 
the Lorentz number attains a universal value L(T) = 
TT 2 kg/3e 2 = Lq. This statement is the essence of the 
Wiedemann-Franz law-i£, formulated long time ago for 
normal metals. Deviations of L{T) from this behavior 
are attributed to non-Fermi liquid behavior—. 

The process used in deriving Eq. ([1]) adopts the some 
of the methods that have been previously developed in 
investigations of the ground state and the transport prop- 
erties of the 2CK modelii^. 



II. THERMAL TRANSPORT 

The minimal Hamiltonian that describes the system 
discussed above is written as H = Hq + iJ int J£. Here, 



r={L,R},a 



dee 4 (e) c r , a (e) , 
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represents the electrons in the non-interacting leads, 
where a single particle state of energy e, spin er and lead 
index r = {R, L} for left, right respectively, is associ- 
ated with the Fermionic operator c\. a (e). The interact- 
ing Hamiltonian contains two terms, H mt = +H^, 
which describe the coupling of the local spin in the 
smaller dot with the conduction electrons in the exter- 
nal leads and to those in the second dot. The first term, 
H {1) is 



H int - \ Jl E 
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Here <x stands for the Pauli matrices er = {<r x , a y , a z } 
and ip^'s are the creation field operators, constructed 
from the creation operators of the electronic states in 



the leads, Vjer = \/6j-D^ ec rcr( e )- Here g is the con- 
duction band density of states, which is electron-hole 
symmetric, g(ui) = 1/2D, —D < ui < D. The cou- 
pling of the local spin with the electrons in the leads 
is considered to be anisotropic, and will be expressed 
through some dimensionless hybridization parameters, 
Vl/r, of the form, r) L / R = v L/R /(v 2 L + v%), with v L/R 
the amplitude of the hopping between the dot and the 
corresponding external lead. An asymmetry parame- 
ter, <j), allows a more elegant description of the hy- 
bridization, til = cos and t/r = sin^^. An 
even/odd basis emerges, with new annihilation operators, 
defined as, * = cos (jfj ip L + sin ( |J ip R and 

^ = sin ( § ) ipL — cos ( f ) ipR. By this unitary transfor- 



mation, the local spin remains coupled only to the even 
channel, while the odd channel becomes decoupled, and 
then, can be treated as a non-interacting one. 
The coupling with the larger dot, 



H int = 9 J 2 E S ^la^oa^la' 



(4) 



is considered to be isotropic, with a dimensionless ampli- 
tude J2. The field operators, ty\ a , describe the electron- 
hole excitations in the larger dot. 

Our theory permits the calculation of the thermal con- 
ductance of a 2CK system by starting from the heat op- 
erator, q(Q>, 

Q iQ) = E / dee(c[je)c L Ae) ~ (4 )(T (e)(c fl , CT (e)) , 

(5) 

which describes the heat transfer under a temperature 
gradient between the leads. We implicitly assume that 
the leads are in equilibrium, fi^ = fiR = 0. The 
heat current is defined through the usual expression 



j(Q) = 



We introduce the so-called com- 



posite fermion operators ; 17 ' 18 T\ = , J ,1 ,er ./ [T S, and 

h) a = 2g 3 / 2 \/3 J^ D deec^(e), to rewrite the heat current 
as 
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The prefactor in the definition of h a was fixed, such 
that it satisfies the canonical anticommutation relations. 
Moreover, h a describes electrons in the odd-channel. 
Since this channel is decoupled it may be treated as a 
non-interacting Fermi system. 

When the system approaches the equilibrium, the heat 
current through the dot can be calculated with the Kubo 
formalism. For that purpose, we consider a temperature 
gradient being applied between the external leads. At 
temperature T, the Hamiltonian acquires an additional 
term, H T = Q (Q) VT/T, with T the temperature itself, 
and VT the temperature gradient between the external 
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leads. Then the average heat current, in its most general 
form, can be expressed as 



(/«>(<))= f Tn{T,t-t')^{t')dt', 

J — OO J 



(7) 



with k[T, t) denoting the thermal conductance. While 
this formalism is general enough to allow the calcula- 
tion of the ac-thcrmal conductance, k(T,u>), here, we 
will focus on the dc-limit only. Algebraic manipulation 
of Eq.([7]) permits expressing the thermal conductance in 
terms of heat current operators, 



Tk(T) 



dt 
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Finally a simple, analytical form can be derived, 



Tk(T) 



2 sin 
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duj u 2 Sm {T{oo. 



T)} 



doj 



(9) 

In Eq. (j9|) we can identify the rescaled Green's func- 
tion of the composite operator as the T-matrix, T(lj) = 

Eq. (jU) allows us to estimate the thermal conductivity 
of the 2CK model. Wc define the Kondo temperature Tk 
through the equation, 3raT(w = Tk) = |3mT(w = 0), 
at the symmetric point, J\ = J2. In this situation, the 
screening of the local spin in the smaller dot is equally 
performed by both channels, their competition leading 
to an NFL quantum critical state. Away from the sym- 
metric point, J\ 7^ J2, the screening is realized by the 
stronger coupling channel. 

In the following, we will always focus on the temper- 
ature dependence of the thermal conductivity for chan- 
nel 1. In the FL regime, for r « T*, the T-matrix is 
analytical^, 



SraTfw, T) ~ — 
irg 



0{K) - sgn(X) 



3u 2 



T' 
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At the symmetrical point, and finite temperatures, the 
non- Fermi liquid character is manifested by a \ T contri- 
bution, 



3m7>, T) ~ (l - ay/u/T K - b^T/T K + ...). 

(11) 

Here a and b are two universal constants. In this case 
the T-matrix is purely imaginary and is reduced by half 
from the unitary value in the Fermi liquid regime. Simple 
integrations leads to the qualitative behavior presented 
in Fig. OD 

Although this is a perturbative calculation, in the FL 
regime, the ratio of the thermal conductance and the elec- 
trical conductance recovers the Lorentz number exactly 
in the T — > limit. It is interesting to point out that 
this also happens for J\ = J2, when again, for liniT->o> 
L(T) = Lq. This result is somewhat contradictory to 
the general conclusion that attributes deviations of the 
Lorentz number from its constant value L to non-Fermi 
liquid behavior. 
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FIG. 2: (color online) Thermal conductance as function of 
T /Tk on a logarithmic scale. The NFL fixed point result, 
Ji = J2, corresponds to the solid line, while dashed and dash- 
dotted lines correspond to Ji / J2. K is the anisotropy 
factor (Ji — J2)/,/ 2 . In both cases the average coupling is 
J = (Ji + J 2 )/2 = 0.2. 
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FIG. 3: (color online) The thermal conductance in the Fermi 
liquid regime for K > (upper panel) and for K < (lower 
panel). k(T)/T given by Eq. |T} is plotted as a function of 
T/T* on a scaled temperature axis. 



III. NUMERICAL RESULTS 

The analytical results obtained in Eq. ([1]) above were 
subjected to a numerical test produced with the numer- 
ical renormalization group (NRG) method. Within the 
NRG framework, the spectral function of the T-matrix 
is given as a weighted sum of (5-functions of the form 
3toT(w, T) ~ Wi5(uj — u>i). This expression replaces 
the integral over frequency in Eq. The numerical 
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FIG. 4: (color online) The universal behavior of the Lorentz 
number as function of T/T* in the Fermi liquid regime. In 
the low temperature limit T <C T* , L(T) converges to Lo. 
The inset presents L(T) as function of T/T K . At the NFL 
fixed point, for T <C Tk the same universality is recovered. 



results were obtained with the Flexible-DMNRG code 1 ^, 
which allows the explicit use of the symmetries of the 
Hamiltonian. In our particular case we built in the 
SU cl {2) eg) SU c2 {2) <g> SU S {2) with SU ca (2) the charge 
SU (2) symmetry in channel a and SU S (2) the global spin 
symmetry. For every run we have kept 1000 multiplets 
and we fixed A = 2. In Fig. [5] we represent the results 
for the thermal conductance k(T) as function of T/Tk- 
At the symmetric point, J\ = J%, the curves for differ- 
ent J's are scaling one on top of each other and we have 
found that the thermal conductance is a universal func- 
tion of T/Tk-, with the behavior predicted by Eq. (JTJ. 
In Fig. [3] we have represented k(T)/T as function of 
T/T*. For temperatures T -C T*, the universal behavior 



k(T)/T cx (T/T*) 2 , obtained analytically in Eq. <JT) was 
also recovered. Thus the numerical results confirm the 
universal behavior found analytically. 

Finally, the universality of the Lorentz number was 
checked numerically. In Fig. 01 main panel, we present 
the results for the Lorentz number as function of temper- 
ature in the FL regime. The conductance G(T) employed 
in the calculation was previously showni^ to have a uni- 
versal behavior for T < T*, when G(T) oc G(T/T*). 
Since both the electric and thermal conductance have 
a similar behavior, L(T) cx L(T/T*) is found to be 
a universal function. At the NFL fixed point, when 

G(T) ~ G (l - V 71 " T/Tk) A th c Lorentz number be- 
comes a universal function of T/Tk- As can be seen in the 
inset of Fig. [U this universal behavior is valid well below 
Tk, while the scaling in the FL regime goes up to temper- 
atures as large as T* . In the limit T — > 0, L(T)/Lq 1 
in all cases. Similarly, a large violation of the WF law 
is observed at large temperatures, T 3> Tk, mostly on 
account of the suppression of the thermal transport. 

In conclusion we have constructed a general frame- 
work for the calculation of the thermal conductance. The 
method can be easily extended to the computation of the 
thermopower, as well as the ac-thermal properties of any 
quantum dot system. The NRG method was used to nu- 
merically evaluate the conductance, thermal conductance 
and the Lorenz number in the 2CK model as function of 
temperature. Our findings point towards a universal be- 
havior of all these quantities, similar with the electrical 
conductance, studied previously^. In the limit of zero 
temperature we have found that L(T — > 0) = Lq both in 
the FL and NFL ground states. 

This research has been supported by Hungarian grants 
OTKA No. K73361, Romanian grant CNCSIS PN II ID- 
672/2008, the EU-NKTH GEOMDISS project, and DOE 
grant number DE-FG02-04ER46139. 



1 D. Goldhaber-Gordon, H. Shtrikman, D. Mahalu, 
D. Abusch-Magder, U. Meirav, and M. A. Kastner, Na- 
ture (London) 391, 156 (1998). 

2 Y. Oreg and D. Goldhaber-Gordon, Phys. Rev. Lett. 90, 
136602 (2003). 

3 T. A. Costi and V. Zlatic, Phys. Rev. B 81, 235127 (2010). 

4 B. Dong and X. L. Lei, J. Phys.: Condens. Matter 14, 
11747 (2002). 

5 J. Koch, F. von Oppen, Y. Oreg, and E. Sela, Phys. Rev. 
B 70, 195107 (2004). 

6 M. Krawiec and K. I. Wysokinski, Phys. Rev. B 75, 155330 
(2007). 

7 A. V. Andreev and K. A. Matveev, Phys. Rev. Lett. 86, 
280 (2001). 

8 K. A. Matveev and A. V. Andreev, Phys. Rev. B 66, 
045301 (2002). 

9 R. Scheibner, H. Buhmann, D. Reuter, M. N. Kiselev, and 
L. W. Molenkamp, Phys. Rev. Lett. 95, 176602 (2005). 

R. M. Potok, I. G. Rau, H. Shtrikman, Y. Oreg, and 



D. Goldhaber-Gordon, Nature 446, 167 (2006). 

11 A. Roman, C. P. Moca, and D. C. Marinescu, unpublised 
(2010). 

12 R. Franz and G. Wiedemann, Ann. Phys. (Berlin) 165, 
497 (1853). 

13 A. Garg, D. Rasch, E. Shimshoni, and A. Rosch, Phys. 
Rev. Lett. 103, 096402 (2009). 

14 M. Pustilnik, L. Borda, L. Glazman, and J. von Delft, 
Phys. Rev. B 69, 115316 (2004). 

15 A. I. Toth, L. Borda, J. von Delft, and G. Zarand, Phys. 
Rev. B 76, 155318 (2007). 

16 I. Affleck and A. W. W. Ludwig, Phys. Rev. B 48, 7297 
(1993). 

17 T. A. Costi, Phys. Rev. Lett. 85, 1504 (2000). 

18 G. Zarand, L. Borda, J. von Delft, and N. Andrei, Phys. 
Rev. Lett. 93, 107204 (2004). 

19 We have used the open access Flexible-DMNRG code, 
available at |http://www.phy.bme.hu/~| dmnrg. 



